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Summary: The purpose of this paper consists in proposing a generalized 
solution for a porous media type equation on a half-line with Neumann 
boundary condition and prove a probabilistic representation of this solution 
in terms of an associated microscopic diffusion. The main idea is to construct 
a stochastic differential equation with reflection which has a solution in law 
and whose marginal law densities provide the unique solution of the porous 
""*■* . media type equation. 
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1 Introduction and preliminaries 

In this work we focus on a porous media type equation given by 

d t u(t, x) € \&l x P (it) (t, x) : (t, x) G (0, T] x R + , 
u(0,x) = u (x) :i£l + , (1.1) 

t 8 X (f3(u))(t, 0) = 0, :te(0,T]. 

The natural analytical concept of (weak) solution of (|1.1|) is given in Def- 
inition H] and it involves the restriction of the derivative of (3(u) on the 
boundary. We introduce here a new notion of solution that we call gen- 
eralized solution for (|1.1|) . which do not require, a priori, the existence of 
distributional derivatives for /3(u). 

Under some minimal conditions, we will first concentrate on uniqueness of 
the generalized solutions (in a large class) and existence of a weak solution 
(smaller class). In particular, we will include the case when /? is possibly 
discontinuous. Moreover we are interested in its probabilistic representation 
through the marginal laws of a stochastic process. 

We formulate now some assumptions. 

Assumption 1 

i) no € L 1 (M+) fl L°° (M+) is an initial probability density; 

ii) j3 : M+ — > M+ is a monotone increasing function with f3 (0) = 0. 

iii) There is a constant c such that \/3 (u)\ < cu. 



With (3 we naturally associate a maximal monotone graph still denoted by 
the same letter /3 : M + — > 2 R +, by filling the gaps, i.e. by identifying (3{x) 
with the interval \p(x—),0(x+)]. We consider now $ : M^_ — > M + such that 
P{u) = & 2 (u)u, u^O. Again we associate naturally $ with the non-negative 



graph (still denoted by $) u i— > y ^p-- Finally we extend the graph $ it 
to R, defining $(0) := [lim inf :r _ > o+ ^ > (x),limsup :r _ !h0+ $(x)]. Of course, if 



(3 is continuous, the first line of (|1-1|) can be replaced by the most natural 
equality dtu(t,x) = g^-L/? (u) (t,x). 

Definition 1. i) We say that (3 is non- degenerate if there is a constant 
Co > such that <£ > Co, /or ewery x > 0. 

m,) We say that (3 is degenerate if lim $ (m) = 0. 

mj M^e say that f3 is strictly increasing after some zero u c if there is u c > 
such that (3 \[q Uc [ = and f3 \\ Uc00 \ is strictly increasing. 

Often also the Assumption 2 below will be in force. 

Assumption 2. We suppose that one of the following properties is verified. 

i) j3 is non-degenerate. 

ii) f3 is degenerate and there is a discrete ordered set of elements ({e^} of 
R+ U {+oo}) so that uq has locally bounded variation on )efc,efc+i(. 

iii) (3 is strictly increasing after some zero u c (in particular it is degenerate) . 

As we anticipated, the idea is to construct a stochastic process Y such that 
the (marginal) law of Yt has a density given by u(t, ■) for any t £ [0, T\. We 
look for Y as being a solution (in law) of the stochastic differential equation 
with reflection 

' dX t = $(u(t,X t ))dB t + dK t , 
K increasing process L Y s dK s = 0, (1-2) 

u (s, .) = law density of Y t , 

which has a weak solution X whose law density is the unique solution of the 
dLH). 

As far as our knowledge is concerned, this paper could be the first one study- 
ing the probabilistic representation of a non-linear partial differential equa- 
tion (PDE) with Neumann boundary conditions on some domain, through 
a reflected non-linear diffusion. 



The problem of probabilistic representation related to a PDE, related to 
a class A of solutions, is the following. For each of u £ A, there exists 
a stochastic process, solving some form of stochastic differential equation 
whose coefficients involve the law of the process (McKean-Vlasov type), 
whose marginal laws are given by u. Solutions of those stochastic differ- 
ential equations are also called non- linear diffusions. The PDE is intended 
as a non-linear forward Kolmogorov's equation corresponding to the non- 
linear diffusion. 

The paper provides a bridge between two big areas of stochastic analysis: 
stochastic differential equations with reflection on some domain, non-linear 
diffusions on the whole line. As we will see the literature is rich of contribu- 
tions in both topics, but in principle no one connects them. 

1. Non-linear diffusion problems 

There are several contributions to the study of equations stated in the 
first line of (|Lip but on the whole line or even on W 1 . That equation, 
which will be precisely stated in (|3.7p . was first investigated by [9] for 
existence, [13] for uniqueness and [10] for continuous dependence on 
coefficients. 

The physical interpretation of the probabilistic representation is the 
following. The singular non-linear diffusion equation (I3.7D describes 
a macroscopic phenomenon for which the probabilistic representation 
tries to give a microscopic probabilistic interpretation via a non-linear 
stochastic diffusion equation modeling the evolution of a single point 
on the layer. 

To our knowledge, the first author who considered a probabilistic rep- 
resentation for the solution of a non-linear deterministic PDE (on the 
whole line), was McKean [22], especially in relation to the so called 
propagation of chaos. He supposed to have smooth coefficients. Af- 
ter that, the literature grew and nowadays there is a vast amount of 
contributions to the subject, particularly when the non-linearity (3(u) 
appears inside the first order part, as e.g. in Burgers equations (see 
for instance the survey papers [H] and [28]). 
A probabilistic interpretation of (|3.7p when j3{u) = \u\u m , for m > 1 



was provided for instance in [Sj in which probabilistic representations 
of the Barenblatt solutions and of a large class of solutions were given. 
Later, when <I> is of class C 3 , Lipschitz, j3 is non-degenerate and uq is 
smooth enough, [19j provided also strong solutions to the probabilistic 
representation problem, see more precisely Remark 1171 In particular, 
the probabilistic representation of the porous media on M was stud- 
ied in the case of irregular coefficients in [11|, 0] with refinements in 
[6] . In particular [11] represented when (3 is non-degenerated all 
the solutions in the sense of distributions under Assumption 1 . More- 
over also the uniqueness of the corresponding non-linear diffusions was 
established. When (3 is degenerate, under Assumptions 1. and 2. ii) 
or 2. hi), [3] has provided again probabilistic representations, but not 
uniqueness of those. Some improvements also appeared in [7, 6], at 
least when <J> is continuous. [7] provides probabilistic representation of 
the Barenblatt solutions when /3(u) = u m , I < m < 1, i.e. in the case 
of fast diffusions. 

2. Stochastic differential equations with reflection. There is a vast 
literature in the subject, in the one-dimensional case, and in the mul- 
tidimensional case as well. It is for us impossible to quote all those. In 
the half- line case M+, such an equation can be formulated as follows: 

dX t = a(t, X t )dW t + b(t, X t )dt + dK t , 

where a, b : [0, T] x 1R are Borel functions and K is an increasing process 
such that J" X s dK s = 0. The solution is the couple (X, K). When 
a and b are Lipschitz the theory is well-developed, see for instance 
|21j . |15j and [25] (and references therein). All the three papers have 
treated the case when the coefficients are time-independent. Elements 
related to the time-dependent case appear for instance in [16] and in 
[23], chapter IX, Exercise (2.14), page 385. 

When the coefficients are non-Lipschitz, but are Lipschitz with a log- 
arithmic correction and the domain is a half- line |12j . proved recently 
the existence of a strong solution. 

As far as weak (in law) solutions of SDEs with reflection, the pioneer- 
ing work is [27], which solves some submartingale problem related to 



stochastic differential equations with reflection under relatively gen- 
eral conditions on the (even time-dependent) coefficients. An interest- 
ing work is also [23J which constructed solutions of symmetric (time- 
homogeneous) diffusions with reflection on a domain of M n , via Dirich- 
let forms. 

Reflected diffusions are naturally candidates for the probabilistic rep- 
resentation of solutions of (linear) PDE with Neumann boundary con- 
dition in the following sense. Given a specific solution u to a Fokker- 
Planck type PDE, with Neumann boundary conditions, a process X 
represents it probabilistically if the marginal laws of X are solutions 
of the PDE. As far as we know, even in this linear case, the point is 
not clear in the literature. For instance, by applying some Ito formula 
type, it is possible to show that solutions of SDEs (or martingale prob- 
lems) solve a PDE with boundary condition, in some sense. However to 
show that a given solution to a Fokker-Planck PDE can be represented 
through a process is rarely explained. This is related to the study of 
uniqueness of the mentioned PDE, see for instance |llj . Theorem 3.8, 
for an equation in the whole line. 

We conclude this discussion about probabilistic representation of lin- 
ear PDEs mentioning [T7| , which gives a representation of the solution 
of an elliptic problem with unbounded monotone drift in term of the 
invariant measure of a reflection diffusion equation considered by [15]. 

We are aware that the Neumann problem on the half-line constitutes some- 
how a toy model, however at our knowledge, there are no results in the 
literature about well-statement of that problem. If x varies in a bounded 
domain (for instance a compact interval), there are some contributions at 
least in the case when /? (u) = u m ,0 < m < 1, which constitutes the case 
of the classical porous media equation. Given an integrable initial condition 
uq, p], through the techniques of maximal accretivity, see e.g. [26], [3], a 
C°-type solution (or mild solution), see Chapter IV. 8 of [2B]. The technique 
consists in showing that the elliptic corresponding operator is m-accretive, 
whose step was performed by [2]. In Corollary 3.5, pQ shows that (when 
uq £ L°°), that the solution is even classical and it is therefore a weak so- 



lution in the sense of Definition HI adapted to the case when an interval 
replaces the real line. 

The paper is organized as follows. After the introduction above, at Section [2] 
we introduce the basic definition of solutions and some notations. At Section 
[31 we discuss existence and uniqueness of (|l-lj) and we remark that the 
solutions have some minimal regularity properties. Section [J] is devoted to 
the existence of the probabilistic representation in the form of a solution to a 
non-linear (in the sense that the marginal densities appear in the coefficients) 
stochastic differential equation, with reflection. The Appendix is devoted 
to the equivalence of weak and generalized solutions under some minimal 
regularity conditions. 



2 Preliminaries 



Let I be a real interval. Given a function ip : [0, T] x i" — > R, (£, x) \- > 
<p(t,x), we denote (if it exists), by cp' (resp. ip") the partial derivative d x <p 
(resp. second partial derivative d xx <p) with respect to the second argument, 
defined again on [0, T] x I. If I is closed, then the derivatives are defined as 
continuous extensions from [0, T] x Intl. 

In this paper C^° (R+) will denote the set of functions <p : M+ — > R which 
are restrictions of smooth functions with compact support defined on R. We 
denote by W^ c (R+) the space of absolutely continuous functions / : M + — > R 
such that for every compact subset K of R + f K \f'(x)\dx is finite. Of course 
for any / G Wj ' c (R+) and every compact K of R we have f K \f(x)\dx < oo. 

If / = R or R+, we denote a bit abusively by WJ£([0,T\ x I) the set of 
/ : [0,T] x R -> I such that for almost all t, we have f(t, •) G W^£(I) and 
for fi 0T i xK \d x f(t,x)\ < oo, for every compact K of /. 

Definition 2. A function u € L 1 ([0,T] x R + ) is called generalized solu- 
tion for equation il-l\) if, for any ip € C^° (R+) such that ip' (0) = 0, we 



have 



t />oo 



ip (x) u (£, x)dx = I tp (x) uq (x) dx + — I I ip" (x) rj u (s, x) dxds, 
Jo ^ Jo Jo 

(2.3) 

where i] u : [0,T] xR + ^f,jj B e L 1 ([0,T] x M+) is such that 

7] u (t, x) G P(u (t, x)) , dt<g>dx - a.e. (t, x) e [0, T] x R + . 

We remark a generalized solution in the same spirit, for porous media equa- 
tions but with Dirichlet boundary conditions was given in [[3] page 226]). 
Moreover, in our case the test function is time independent. 

Remark 3. The formal justification of previous formula comes out from the 
following observation. Suppose that (u, rj u ) is a smooth solution to II 1.1]) and 
uq is continuous. Then, for every x G M +; we have 

1 /"* 

u(t,x) = u (x) + - r]"(s,x)ds. 
z Jo 

Let ip be a test function from C^° (M+) such that ip' (0) = 0. In this particular 
situation we have 

OO ^00 -1 pt pOQ 

ip(x)u(t,x) dx = / tp (x) Uq (x) dx H — / / <p (x) ry" (s , x) dxds 
o Jo ^ Jo Jo 

OO 1 ft 

ip (x) u (x) dx + - / <p (0) ^ (s, 0) ds 
ii ^ Jo 

£ /"OO 



- / / y/ (x) rf u (a, x) dxds 

1 Jo Jo 



OO 1 ft 



p (x) uq (x) dx + - ip (0) 77^ (s, 0) ds 



I! 2 j 



1 /■* f 00 1 /■* 

+- / / ip" (x)r] u (s,x)dxds- - <p (0)r] u (s,Q)ds, 

1 Jo Jo l Jo 

y v ' 

=0 

which constitutes indeed (|2.3p . 

As we mentioned in the introduction, the natural analytical concept of so- 
lution should involve the first spatial derivative at the boundary. 



We consider a couple (u,rj u ) satisfying Definition [2] and suppose that r\ u G 
Wj Q ' c ([0, T] x R_|_) in agreement with the notations of Section [2j In particular 
for every compact K £ R we have 

/ \r)' u (s,x)\ dsdx < oo. (2.4) 

J[0,T]xK 

Definition 4. The couple (u,rj u ) is said to be a weak solution of (|l.lj) if for 
every tp € Co°(R + ) we have 

oo poo pt /"OO 

(p (x) u (t, x) dx = (f (x) uq (x) dx — ip' (x) r\ u (s, x) dxds. 

o Jo Jo Jo 

(2.5) 

Proposition 5. Let (u,rj u ) be a couple such that r/ u (t,x) € f3(u(t,x)) and 
f] u € W^ ' c ([0, T] x M+). Then u is a generalized solution if and only if it is 
a weak solution in the sense of Definition^ 

Proof. See Appendix. □ 



3 The porous media equation on half-line with Neu- 
mann boundary condition 



In this part of the paper we will study existence and uniqueness of the 
generalized solution for equation (jl.ip . We also show the connection with 
the notion of weak solutions. 

This will be done using the known results on the whole line M.. To this 
purpose, we start extending the initial condition to the real line by the 
following construction. 
Let Ho- £ (L 1 n L°°) (R) be defined by 

«oW= / , . n (3-6) 

[ 2 M ° \ X ) , 37 < 0. 

and p : R -> R by 

P(u) = -/3(2u), ueR. 



We can now consider the corresponding porous media equation on the whole 

line 

d t u(t, x) = \dlJ5 (u) (t, x) : (t, x) G (0, T] x R, 

(3.7) 
u(0, .) = u , 

which, by Proposition 3.4 from [II] (see also [10]) has a unique solution 
in the sense of distributions, i.e. there exists a unique couple (u, rj u ) G 
(L 1 n L°°) ([0,T] x R) such that 

tp (x) u (t,x) dx = / ip (x) uq (x) dx+- / / (p" (x)r]u(s,x) dxds, (3.8) 
for all <p G C°° (M) and 

r]u (t, x) G /3 (u (£, x)) , for d£ <g> dx - a.e. (t, x) G [0, T] x R. 



Remark 6. Since Uq is even, the solution u to equation \3. 7| ) and ra^ are ako 
even. In fact, by applying (|3.8j) to x h- > </? (— x) and making then the change 
of variables x i-> — x, we s/iou> £/ia£ (£,x) i-> u(£, — x) is a/so a solution. The 
result follows by uniqueness of (|3.7p . 



Proposition 7. VFe define v,r] v : [0, T] x R +; setting v(t,x) = 2u(t,x) and 
r] v (t,x) = 2r]u(t,x),\/(t,x) G [0,T] x R + . T/ie couple (v,r] v ) is a generalized 
solution to equation hl.l\) in the sense of Definition^ 

Proof. Since 

we first observe that rj v (t,x) G (3{v{t,x))dtdx a.e. Moreover rj v belongs to 
L 1 ([0, T] x R + ) since r/ s belongs to L 1 ([0,T] x R). It remains to prove that 
(v,rh) satisfies (03]) for all ip G C^ (R+) such that <// (0) = 0. 

Step I 

First we prove that (J2.3J) is true for all test functions </> G C^° (R+) with 
<// (0) = and which extend to an even smooth function (p with compact 
support on R. By (j3.8|) we have 



/ 

Ju 



2./o 



93 (x) u (t, x)dx -- I ip (x) uq (x) dx + - / / ¥>" (&) ry« (s, x) dxds. 
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Since p, uq and u(t, •), for any t G [0, T] are even functions, we get, for the 
first two terms, that 

/ Tp(x)u(t,x) dx = / tp(x)2u(t,x)dx = / tp (x) v (t, x) dx, 

Jr Jo Jo 

/ Tp (x)uq (x) dx = j Tp (x) 2uq (x) dx = / Tp (x) uq (x) dx . 

Jr Jo Jo 

Since rj v = 2r] u we obtain 



oo 



ip" (x) rju (s, x) dx — j ip" (x) 2r\u (s, x) dx 

tp" (x) rj v (s, x) dx, s G [0, T] 



o 

oo 



/() 

This proves (jl.lj) for the restricted class of </? which extend to an even smooth 
function on R with compact support. 

Step II 

Now we can prove the general statement. 

Let ip € C^° (R+) such that <// (0) = 0. We extend it to an even function 

_ f (p(x) :x>0 

[ ip{— x) : x < 0, 

which has compact support, but it does not belong necessarily to Cq° (R) . 
In order to have a proper test function for evaluating it in (|3.8p , we need to 
convolute it with a mollifier. 

Let p G C^° (R) be such that p > for |x| > 1, p(x) = p(— x) and 
L p (x) dx = 1. For an example of such function see e.g. [2]. 

We set /3 £ (x) = —p I — j as a mollifier and we take the regularization 

Ve (x) = hp (x) pe {X ~ y) dy 

Jr 

Tp(x — ey) p (y) dy, Vx G R. 

It is well-known that Tp e G C^° (M) and we can check that Tp e is also even. 



11 



Indeed we have 



<p e (-x) = / y{-x-ey)p{y)dy 
Tp (x + ey) p (y) dy 
Tp(x- ey) p (-y) dy = Tp £ (x) 



By Step I we get 

fOO 1 ft /'OO 

Tp e (x) v (t, x) dx = j Tp £ (x) uq (x) dx+- / / (jp e )" {x)rj v (s,x) dxds. 
Jo % Jo Jo 

Since Tp £ L + G C^° (R+) we have 

&e)" ( x ) = / ¥' ( x - £ y) p (y) d y = (jp") e ( x ) > Vi g r+; 

then we can pass to the limit for e — > and conclude the proof. □ 

Corollary 8. Under Assumptions 1. and 2., the generalized solution of 
(II. ip is also a weak solution. 

Proof. Let (u,r]u) be the solution of (|3.7p . i.e. on the real line with initial 
condition tio as i n (|3.6p . Under Assumptions 1. and 2., Proposition 4.5 a) of 
[4] says that for a.e. t G [0,T], r/ n (£, •) G iT^R) and /r rixM r ?«( i ' x ) 2 dtdx < 
oo. This implies that for every compact real interval X, ry^ G L 1 ([0,T] x iiT). 
In particular r/u G W^ ' ([0, T] x R). By the proof of Proposition [9j the 
solution (u,T] u ) equals the restriction of (2u,2fj) to R + . This shows that rj u ) 
belongs to W^ ' c ([0, T] xR + ). By Proposition^ u is also a weak solution. □ 

Proposition 9. Equation U.l\) has a unique generalized solution in the 
sense of Definition^ 

Proof. Existence has been the object of Proposition so we proceed now 
to uniqueness. Let (v,r] v ) be a generalized solution of (|1.1|) . i.e. for any 
^eC °° (R+) such that iff (0) = we have 



OO fOO 1 ft fCG 

tp(x)v (t,x) dx = / ip (x) uq (x) dx + - / / 
i) Jo % Jo Jo 



ip (x) v (t, x) dx = j ip (x) uq (x) dx + - / / ip" (x) r} v (s, x) dxds. 

2 Jo Jo 

(3.9) 
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We define u : R -> R by 

«o(x) 



^«o (x) , x > 



i ^«o (— x) , x < 
We extend the considered solution to [0, T] x R setting 

_ j |u(t,x) ,x>0 

u(t,x) = \ , 

y 7jt> (£, — x) ,x < 0, 

,, s / 5% (*> s) , £ > 

% (*, X) = < ! , . , 

[ %q v {t,—x) ,x < 0. 

Obviously we have 

% (*, a) € -/? (2u (£, a;)) = J3(u(t, x)) ; dt®dx- a.e. (t, x) 6 [0, T] x R + . 

We aim now at showing that («,%) constructed above, is a solution to 
equation (|3.7|) in order to use the uniqueness for the equation on the whole 
line. To this purpose we have to prove that 

tp(x)u(t,x) dx = / (p (x) uq (x) dx + - / / if" (x)i]u(s,x) dxds, 
Jul 2 Jo 7r 



for all tp G C^ 

Let us now fix ip € C^° (R) and we can define ip : R+ — > R by tp (x) = 
c/?(x) + (f(-x). Obviously ip G C^ (R + ) and V' (0) = 0, so it is a test 
function for (|3.9|) . This gives 



ip(x)v (t,x)dx = / -0 (x) uo (^) drc + - / / ip" (x)r] v (s,x)dxds. 
IQ Jo 2 ,/ Jo 

The left-hand side above yields 

00 /"OO 

ip (x)v (t,x)dx = / (y? (x) + y? (— x)) v (t, x) dx 

I! Jo 

OD r0 

(f (x) v (t, x) dx + / ip (x) i> (£, — x) dx 

J-oo 



tp (x) 2u (t,x) dx, 
by the obvious change of variable x 1— > —x. The same technique shows that 

r _ 

ip (x) Uq (x) dx = if (x) 2lio (x) dx, 

./k 
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and also 

fOO 



/*oo P 

I ifi" (x) r] v (s, x) dx = 2 / ip" (x) r\u (s, x) dx, for all s G [0, T] . 
Jo Jr. 

This implies 

2 / ip (x) u(t, x) dx = 2 / ip (x) uq (x) dx + 2 / / 99" (x)r]u(s,x) dxds, 

Jr Jr Jo Jr 

where 

7]u (t, x) e J3 (u (t, x)) , a.e. (t, x) £ [0, T] x R + 

and we get that 

J d t u{t,x) € \dl x {fi{u{t,x))) :(t,x)e(0,T]xR, 
\ u(0,x) = u Q -( x ). 



Finally, this shows that u is a solution to equation (I3.7p . As said at the 
beginning of the proof, the uniqueness for v follows from uniqueness of (|3.T|) . 

□ 

Remark 10. By Remark 1.6 of [$, the solution u of (|3.7p on the real line, 
has mass conservation and it is always non-negative. Taking into account 
the construction of the generalized solution, we also get a similar result for 
the generalized solution v of (jl.ip . This means the following. 

1. v > a.e. 

2. f R v(t,x)dx = 1, ViG [0,1]. 

This explains why the values of f3 on M_ are not important, see Assumptions 
1. and 2. 



4 The basic construction of the linear reflected dif- 
fusion 



In this section we are interested in the probabilistic representation of equa- 
tion (jl.ip . under Assumptions 1. and 2. from the introduction. More 
precisely we aim at characterizing all the solutions of (jl.ip . 
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In this case we remind that we can write j3 (u) = 5> 2 (it) u, for $ being a 
non- negative graph generated by a bounded function $. 

We introduce now <J> : M — > R by $ (it) = $ (2u) . With this notation we have 
j3 (u) = u $ (it) . Let mq be as in (|3.6[) . 



According to Theorem 4.4 in [11] , Theorem 5.4 in [1] taking into account 
the proof of Theorem 2.6 in [5], equation (|3.7p admits a probabilistic repre- 
sentation, in the sense that there is a solution Y (in law) of 

' Y t GY + J * $ (II (a, y s )) cW s , t G [0, T], 
H(£, .) = law density of Y t ,t G [0,T], (4.10) 

u(0, •) =« 0) 

for some Wiener process (W s ) on some probability space. 
The precise meaning of the first line of (|4.10p is the following: 
Y t G y + / Xu(s,Ys)dW s ,t£ [0,T], x *(t,v) G $(u{t,y)),dtdy a.e. (4.11) 

JO 

In order not to loose the reader, without restriction of generality, we suppose 
in this section that <3? (and therefore $) is continuous on M_|_ so that we can 
write Xu = $(«)• 

Remark 11. VFe remind that we have defined $ ("and therefore <&) artifi- 
cially at zero in the Introduction. Observe that this has no influence in the 
probabilistic representation since L lt a f s Y s )=o\^ 2 (u(s,Y s ))ds = a.s. since 
its expectation gives J ds f^lc U ( s ^ =0 }u(s,y)Q 2 (u(s,y))dy = 0. 

Our purpose is to use this result in order to get a corresponding one for the 
porous media equation on the half-line with Neumann boundary condition. 
We start with a preliminary result. 

Lemma 12. Let (Y,u) be the solution of fl^. ICfy . Let us denote by vt,t G 
[0, T], the marginal laws of Xt = \Yf\,t G [0, T]. Then for all t G [0, T], ut 
has a density which is given by v(t, •) where v = 2 H|[o,t]xR + - 
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Proof. We remind that u is the law density of Yt from equation (|4.10p . Let 
cp : M+ — > R be a bounded Borel function and t € [0, T\. We notice that 

M(<p(X t )) = E(<p(\Y t \))= [ <p(\y\)u(t,y)dy = 
r0 



/ 



ip(-y)u(t,y)dy+ / ip (y) u (t , y) dy 



oo 
oo 



v (y) 2 « (*, y) ^y; 



since u(t, •) is even, this concludes the proof. 



□ 



Let v be as in Lemma [T2l We continue by investigating the stochastic 
equation solved by Xt = \Yt\ , t € [0, T] where Y is the solution of (|4,10p . 

By Ito-Tanaka formula (see e.g. [24J Theorem (1.2) Chapter VI) we get 

X t = \Y t \ = \Y \ + f sgn (Y s ) $ (u(s,Y s )) dW s + L Y t (0) , (4.12) 

Jo 

where 

1 ,x>0 

sgn{x) = I —1 , x < 

,x = 

and (Lj (0)) is the local time of the semimartingale Y at zero. Since Y is a 
local martingale, this is characterized by 

i ft 

Y 



U (o) 



sSoi i 1 ^^-^)}^^},. ( 4 - 13 ) 



see e.g. Corollary (1.9) Chapter 6 [21]. We define 

?t 



B\ 



sgn {Y s ) dW s 



o 



B 



ids. 



so that 

ft pt 

(sgn) 2 (Y s ) ds = / l{y s?;i j' 
/o Jo 

Possibly enlarging the probability space, let B 2 be an independent Brownian 
motion of W. We set 



B := B 1 + 



o 
16 



l{ Ya =o}dB 2 . 



(4.14) 



Since (W, B 2 ) t = 0, it follows that 

(B) t = f l {Y ^ 0} ds + / l {Ys=0} ds = t,Vte [0,T\. 
Jo Jo 

By Levy characterization theorem of the Brownian motion, it follows that 

B is a Wiener process. 

By (14.12p . since L is a bounded variation process, we show that 

(X) t = [ l {Y3m <5 2 (u( S ,Y s ))ds= [ l {Yam $ 2 (2u(s,Y s ))ds (4.15) 
Jo Jo 

ft 



y 

because u(s,y) = —(s, \y\), s E [0, T], i/£R. Coming back to (|4.13p . again 
by Corollary (1.9) Chapter VI of [24\ and (J4TTU]) . we get 

L " (0) = iTohl! 1 ^-^^ 

1 /•* _ 2 _ 

= J™^ / l{y s e(-e,e)}^ (u(s,y s ))ds 

I ft 

= lim ^/ 1 {|n|e[o,e)}* 2 (2u(«i*i))<k 

£-S>0 

1 /"* 

= lim ^/ 1 {x s e[o, £ )}^ 2 (^(s,^s))^ 

£->0 

1 /** 

= lim ^/ 1 {x s e(o, £ )}^ 2 (w(s,^s))ds 

£->0 

1 /"* 

+ lim— / l {Xa=0} $ 2 Ms,X s ))d S . 
because of Lemla [12j By (|4.15p it follows that 



1 /"* 
+Jim- y l { x a =o}$ 2 (v (s, X,)) ds. 



(4.16) 



Now 



/ l {Xs=0} <S> 2 (v(s,X s ))ds = 0, (4.17) 

Jo 



'0 
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because the expectation of the non-negative left-hand side gives 

/ / <5> 2 (v(s,y))v{s,y)dyds = 0. 
Jo J{0} 

By ()4.17p . (|4.16|) and taking also into account again Corollary (1. 
VI of [M] we obtain that 



x 



L\ (0) = -Lf (0) 



Going back to (|4.12|) we obtain 

' x t = x + ; * $ (v (s, X s )) dBl + hf (0) 

v (t, .) = law density of Xt, t € [0, T] , 
k v(0,.) = u . 

Taking into account (|4.14p we obtain 

X t = X + f* d> (v (s, X t )) dB s + ~Lf (0) ; 
v (t, .) = law density of X t , 
v(0,-) = u . 



Chapter 



(4.18) 



(4.19) 



This happens because J Q $(v(s, X s ))dB 2 = 0. In fact its expectation gives 
E (f <fr 2 (v(s,X s ))ds) = since for each s € [0, T], X s has a density. The 
proof is now concluded. 



5 The probabilistic representation 



A solution of (J4.19J) is the candidate for a probabilistic representation of a 
solution to (|4.18p . The proposition below gives an illustration of this, even 
though this constitutes the easy part of the probabilistic representation. 

We suppose again the validity of Assumptions 1. and 2. 

Proposition 13. Let X be a solution in law of 

' X t = X + Jl X v (s, X„) dB s + h? (0) , 
v (t, .) = law density of Xt, t € [0,T] , 
Xv (s,x) G $ (v (s,x)) , ds®dxa.e. (5.20) 

X >0, 

k v{0,-) =u . 
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Then v is the generalized solution in the sense of Definition to equation 

Remark 14. i) Equation ft 5. 20]) is formalized by 

X t GX + Jj $ („ ( s , X ,)) dB s + hf (0) , 
u (t, .) = Zaw density of X t , t € [0, T] , u(0, •) = 0, ( 5 - 21 ) 

x X >0. 

m^) //<£ is continuous 115.21]) . coincide with \^.19 ). 

Hi) According to Corollary i> is ako i/ie lueafc solution of 1^4-18 ) in the 
sense of Definition^ 



Proof. Let ^C m i 



be such that ipf (0) = 0. By ltd formula we obtain 



<p(X t ) = <p(X )+ ip'(X s ) X v(s,X s )dB s (5.22) 

Jo 

+\ f ip"(X s )xUs,X s )ds+ f <p'(X s )dL?(0). 
z Jo Jo 

The last term above gives 

/ ^'(X s )dLf(0)= f l {X3=0} <p'(X s )dL?(0)+[ l {x , m <f/(X a )dL?(0). 
Jo Jo Jo 

The first integral of the right-hand side equals 

¥>'(0) / l {Xs =o}dLf (0) = 0, 
Jo 

by the assumptions on ip. 

The absolute value of the second integral is bounded by 



ML/ Wo}<&f (o) = o, 

Jo 



by Theorem 7.1 from [20J. So we can conclude that f Q tp' (X s ) dL^ (0) = 0. 
Now, by taking the expectation in (|5.22p . for t € [0, T], we obtain 



tp (x) v (t,x) dx 



(p (x) uq (x) dx 



+-. 



t poo 



JO 



(p" (x) v (s, x) xt (s, x) dxds. 
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Since /3 (u) = u & 2 (u) , setting r/ v (s, x) = v (s, x) Xv ( s > x ) > we have rj v (s, x) € 
j3 (v (s, x)) , ds®dx a.e. and v is a generalized solution. We remind that v is 
necessarily the unique solution of (jl.ip and we can conclude the proof. □ 



Remark 15. Any solution to equation 



dX t £ $ (v (t, X t )) dB t + \dLf (0) , 
v (t, .) = law density of X t , t G [0,T] , 
v(Q, .) = n , 



solves also the reflecting Skorohod problem 



dX t e<f>(v(t,X t ))dB t + dK t , 

K is increasing, 

X > 0, J T X s dK, = 0, 



(5.23) 



wto Lf (0) = 2K t . 

Indeed, {Lf (0)) is an increasing process, being a limit of increasing pro- 
cesses (see Corollary 1.9 Chapter VI of l2$ ) and f l{x s =o}dLlf (0) = by 
Theorem 7.1 from \20§ . 

In particular showing the existence of a solution to (I5.2ip . also establishes 
the existence for equation [57i 



The consideration of the first part of Section 3 allows to state effectively the 
existence of solutions for (|5.20p which consists in a sort of converse statement 
of Proposition Q31 

Theorem 16. Under Assumptions 1 and 2 there is a solution in law of 

Proof. The proof is the object of Section H] with some obvious adaptations 
to the case when $ was not continuous and it has to be associated with a 
graph. 

□ 
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Remark 17. The conclusion of Theorem [76| is valid under other settings 
of similar hypotheses (but different) as those in Assumptions 1 and 2. For 
instance under the one of the following assumptions, there is a solution in 
law to ([OUj) . 

• Suppose <J> to be continuous on M+ and Assumption 1. ii) is valid. 
Then the conclusion of Theorem{T^ still holds if in Assumption 1. we 
replace uq € L 1 (M) n L°°(K) with uq € L 1 (M). In particular, there is 
a solution in law to to (|5.20p . Indeed by Theorem 2.6 of [5], a prob- 
abilistic representation on the whole line takes place. This constitutes 
the replacement tool to Theorem 4-4 °f IHj an d °f Theorem 5.4 of J^]/, 
mentioned in the lines before (|4.10|) . 

• Item Hi) in Assumption 1. is somehow technical and it can be often 
relaxed for instance when (3 is non- degenerate and <3? is smooth. Typi- 
cally, when $ is of of class C 3 , Lipschitz, (3 is non- degenerate and uq 
is absolutely continuous with derivative in H 2+a for some < a < 1. 
Then there is a (even a strong) solution of (J4.10J) . see for instance 
Proposition 2.2 in \lty - As for the previous item, by the same proof as 
for Theorem \lb\ there will be a solution to (|5.2U|) . 

6 Appendix 

Proof of Proposition [5]. 

Let ip & C^°(IR + ). Operating by integration by parts, (|2,5p it is equivalent 
to 

OD POD 

(f (x) u (t, x) dx = j (p(x)uo{x)dx (6.24) 

o Jo 

t ft i-oo 

if' (0) r?u (s, Q)ds+ / ip" (x) r] u (s, x) dxds. 
o Jo Jo 

We observe that (u,r] u ) is a generalized solution if (|6.24p holds for every 
<p G C^° (M + ) such that iff (0) = 0. 

We can easily see that if (u,r] u ) is a weak solution then it is obviously also 
a generalized solution. 
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In the other sense we assume that (u, %) be a generalized solution. To show 
that it is also a weak solution we need to show that (|2.5p holds for every 
^ G Cfi° (R + ) . We know in fact that (ET22J) holds for every ^ e C °° (R + ) 
such that if' (0) = 0. 

Let e > 0. We consider a sequence Xe 0*0 '■ R+ - ► [0, 1] smooth such that 

( s Jo :xe[0,£) 
1 : x > 2e. 



We set 



where 



Ve (a?) = Xe (y) <ft (y) dy + <p(0) + c (e) 



ii 



c(e) = / (1 - Xe) (l/) ^ (y) dl/> 

JO 

We note that <^ e is obviously smooth. Moreover it has compact support since 
for x > max(supp(^) we have 

/•oo 

Ve (x) = / Xe (y) f' (y) dy + ¥3(0) + c (e) = 0. 

J0 

By assumption, (|6.24|) holds with (^ replaced by ip e . Obviously ip e — > (p 
pointwise and <p' e — >■ <p' in LP (M) , for all p > 1. 

Since <p' e (0) = 0, f)2.5[) holds for y? replaced by ip £ , i.e. 

oo /-oo 

(p £ (x) u (t, x) dx = / ip e (x) uq (x) dx (6.25) 

n Jo 

t POD 

/ ip' e (x)r]' u (s,x)dxds. 
o Jo 

We remark that c{e) — > for e — > and 

/•OO 

|Ve(x)|< / |(/p'(y)|dy + (/j(0) + sup c(e). 

J0 e 

Let i € [0)^1- By the dominated convergence theorem the left-hand side 
of (J6.25P (respectively the first integral in the right-hand side), converges to 
J °° tp (x) u (t, x) dx (respectively / °° ip (x) uq (x) dx). 

We observe 

\<Pe(x)\ < \(ff (x)\li(x), 
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/ being a compact interval including the support of (p. 

Again by dominated convergence theorem the second integral in the right- 
hand side of (16.25f) . converges to J Q j tp' (x) rj' u (s, x) dxds. This shows (12. 5ft 
for (jl.ip and conclude the proof. 
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